Set up
suppressPackageStartupMessages({
library(tidyverse)
})
Read in data and process
pbta_df <- readr::read_tsv(pbta_file, guess_max = 100000, show_col_types = FALSE) %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, cg_multiple, cg_id, cgGFAC, tumor_descriptor)
tmb_vaf_df <- readr::read_tsv(tmb_vaf_file, guess_max = 100000, show_col_types = FALSE) %>%
filter(!tmb >= 10) %>%
select(Kids_First_Biospecimen_ID, Variant_Classification, gene_protein, mutation_count, region_size, tmb, VAF)
genomic_paired_df <- readr::read_tsv(genomic_paired_file, guess_max = 100000, show_col_types = FALSE) %>%
left_join(pbta_df, by = c("Kids_First_Participant_ID")) %>%
left_join(tmb_vaf_df, by = c("Kids_First_Biospecimen_ID")) %>%
filter(!is.na(tmb))
# Attention as some bs specimen might not have TMB!
# If that happens, we will end up with samples lacking timepoints.
# Which patient samples don't have TMB?
# genomic_paired_df %>%
# filter(is.na(tmb)) %>%
# unique() %>%
# regulartable() %>%
# fontsize(size = 12, part = "all")
descriptors_df <- genomic_paired_df %>%
group_by(Kids_First_Participant_ID) %>%
summarize(descriptors = paste(sort(tumor_descriptor), collapse = ", "),)
# Vector to order timepoints
timepoints <- c("Diagnosis", "Progressive", "Recurrence", "Deceased", "Second Malignancy", "Unavailable")
# Create df to use for plots
df_plot <- genomic_paired_df %>%
left_join(descriptors_df, by = c("Kids_First_Participant_ID", "descriptors")) %>%
mutate(td_cgGFAC = case_when(grepl("Deceased", tumor_descriptor) ~ "xDeceased",
TRUE ~ tumor_descriptor),
match_id = paste(tumor_descriptor, Kids_First_Participant_ID, sep = "_"),
tumor_descriptor = factor(tumor_descriptor),
tumor_descriptor = fct_relevel(tumor_descriptor, timepoints))
# Let's count number of samples
count_df <- df_plot %>%
group_by(tumor_descriptor, cg_id, match_id, Variant_Classification) %>%
dplyr::count(cg_id)
Define parameters for plots
# Read color palette
palette_df <- readr::read_tsv(palette_file, guess_max = 100000, show_col_types = FALSE)
# Define and order palette
palette <- palette_df$hex_codes
names(palette) <- palette_df$color_names
Alterations per timepoint
# Define parameters for function
x_value <- count_df$tumor_descriptor
title <- paste("Variant types in PBTA cohort", sep = " ")
# Run function
fname <- paste0(plots_dir, "/", "Alteration_type_timepoints_barplots.pdf")
print(fname)
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/Alteration_type_timepoints_barplots.pdf"
p <- create_stacked_barplot_variant(count_df = count_df, x = x_value, palette = palette, title = title)
pdf(file = fname, width = 6, height = 6)
print(p)
dev.off()
quartz_off_screen
2

Alterations per timepoint in each cancer type
cg_id_id <- as.character(unique(count_df$cg_id))
cg_id_id <- sort(cg_id_id, decreasing = FALSE)
cg_id_id
[1] "Adamantinomatous Craniopharyngioma" "Atypical Teratoid Rhabdoid Tumor" "Chordoma"
[4] "Choroid plexus carcinoma" "CNS Embryonal tumor" "Craniopharyngioma"
[7] "Diffuse midline glioma" "Dysembryoplastic neuroepithelial tumor" "Embryonal tumor with multilayer rosettes"
[10] "Ependymoma" "Ewing sarcoma" "Ganglioglioma"
[13] "Glial-neuronal tumor" "Hemangioblastoma" "High-grade glioma"
[16] "Low-grade glioma" "Malignant peripheral nerve sheath tumor" "Medulloblastoma"
[19] "Meningioma" "Neuroblastoma" "Neurofibroma/Plexiform"
[22] "Pilocytic astrocytoma" "Rosai-Dorfman disease" "Schwannoma"
# Define parameters for function
x_value <- count_df$match_id
title <- paste("Variant types in PBTA cohort across cancer groups", sep = " ")
# Run function
fname <- paste0(plots_dir, "/", "Alteration_type_timepoints_cg_id_barplots.pdf")
print(fname)
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/Alteration_type_timepoints_cg_id_barplots.pdf"
p <- create_stacked_barplot_variant_cg_id(count_df = count_df, x = x_value, palette = palette, title = title)
pdf(file = fname, width = 30, height = 18)
print(p)
dev.off()
quartz_off_screen
2

Alterations per timepoint in each cancer type defined by cgGFAC
df_plot_cgGFAC <- df_plot %>%
group_by(tumor_descriptor, cgGFAC, Kids_First_Biospecimen_ID, Variant_Classification) %>%
dplyr::count(cgGFAC)
cgGFAC_id <- as.character(unique(df_plot_cgGFAC$cgGFAC))
cgGFAC_id <- sort(cgGFAC_id, decreasing = FALSE)
cgGFAC_id
[1] "ATRT" "DMG" "HGG" "LGG" "Other"
# Define parameters for function
x_value <- df_plot_cgGFAC$match_id
Warning: Unknown or uninitialised column: `match_id`.
title <- paste("Variant types in PBTA cohort across cgGFAC", sep = " ")
# Run function
p <- create_stacked_barplot_variant_cgGFAC(count_df = df_plot_cgGFAC, x = x_value, palette = palette, title = title)
Error in `geom_bar()` at ]8;line = 173:col = 2;file:///Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/util/function-create-barplot.Rtmb-vaf-longitudinal/util/function-create-barplot.R:173:2]8;;:
! Problem while setting up geom.
ℹ Error occurred in the 1st layer.
Caused by error in `compute_geom_1()`:
! `geom_bar()` requires the following missing aesthetics: x
Backtrace:
1. global create_stacked_barplot_variant_cgGFAC(...)
3. ggplot2:::print.ggplot(...)
5. ggplot2:::ggplot_build.ggplot(x)
6. ggplot2:::by_layer(...)
13. ggplot2 (local) f(l = layers[[i]], d = data[[i]])
14. l$compute_geom_1(d)
15. ggplot2 (local) compute_geom_1(..., self = self)

Alterations per timepoint in each cancer type and timepoint
model
tm_df_plot <- df_plot %>%
filter(!is.na(timepoints_models)) %>%
group_by(tumor_descriptor, cg_id, timepoints_models, Kids_First_Biospecimen_ID, Variant_Classification) %>%
dplyr::count(timepoint_cgGFAC_n)
tm <- as.character(unique(tm_df_plot$timepoints_models))
tm <- sort(tm, decreasing = FALSE)
tm
[1] "Dx-Dec" "Dx-Pro" "Dx-Pro-Dec" "Dx-Pro-Rec" "Dx-Pro-Rec-Dec" "Dx-Rec" "Dx-Rec-Dec" "Dx-SM"
[9] "Pro-Dec" "Pro-Rec" "Pro-Rec-Dec" "Rec-Dec" "Rec-SM"
# Loop through variable
for (i in seq_along(tm)){
print(i)
df_sub <- tm_df_plot %>%
filter(timepoints_models == tm[i])
# Define parameters for function
x_value <- df_sub$Kids_First_Biospecimen_ID
title <- paste(tm[i])
# Run function
p <- create_stacked_barplot_variant_cg_id(count_df = df_sub, x = x_value, palette = palette, title = title)
}
[1] 1
[1] 2
[1] 3
[1] 4
[1] 5
[1] 6
[1] 7
[1] 8
[1] 9
[1] 10
[1] 11
[1] 12
[1] 13













LS0tCnRpdGxlOiAiQ2xhc3NpZmljYXRpb24gb2YgVmFyaWFudHMgYWNyb3NzIHBhaXJlZCBsb25naXR1ZGluYWwgc2FtcGxlcyBpbiB0aGUgUEJUQSBDb2hvcnQiCmF1dGhvcjogJ0FudG9uaWEgQ2hyb25pIDxjaHJvbmlhQGNob3AuZWR1PiBmb3IgRDNCJwpkYXRlOiAiMjAyMyIKb3V0cHV0OgogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IFRSVUUKICAgIHRvY19mbG9hdDogVFJVRQotLS0KCiMgU2V0IHVwCmBgYHtyIGxvYWQtbGlicmFyeX0Kc3VwcHJlc3NQYWNrYWdlU3RhcnR1cE1lc3NhZ2VzKHsKICBsaWJyYXJ5KHRpZHl2ZXJzZSkKfSkKYGBgCgojIERpcmVjdG9yaWVzIGFuZCBGaWxlIElucHV0cy9PdXRwdXRzCmBgYHtyIHNldC1kaXItYW5kLWZpbGUtbmFtZXN9CiMgRGV0ZWN0IHRoZSAiLmdpdCIgZm9sZGVyIC0tIHRoaXMgd2lsbCBiZSBpbiB0aGUgcHJvamVjdCByb290IGRpcmVjdG9yeQojIFVzZSB0aGlzIGFzIHRoZSByb290IGRpcmVjdG9yeSB0byBlbnN1cmUgcHJvcGVyIHNvdXJjaW5nIG9mIGZ1bmN0aW9ucyAKIyBubyBtYXR0ZXIgd2hlcmUgdGhpcyBpcyBjYWxsZWQgZnJvbQpyb290X2RpciA8LSBycHJvanJvb3Q6OmZpbmRfcm9vdChycHJvanJvb3Q6Omhhc19kaXIoIi5naXQiKSkKYW5hbHlzaXNfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImFuYWx5c2VzIiwgInRtYi12YWYtbG9uZ2l0dWRpbmFsIikKcmVzdWx0c19kaXIgPC0gZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgInJlc3VsdHMiKQppbnB1dF9kaXIgPC0gZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgImlucHV0IikKZmlsZXNfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImFuYWx5c2VzIiwgInNhbXBsZS1kaXN0cmlidXRpb24tYW5hbHlzaXMiLCAicmVzdWx0cyIpCgojIElucHV0IGZpbGVzCnBidGFfZmlsZSA8LSBmaWxlLnBhdGgoZmlsZXNfZGlyLCAicGJ0YS50c3YiKSAjIGZpbGUgZnJvbSBhZGQtc2FtcGxlLWRpc3RyaWJ1dGlvbiBtb2R1bGUKZ2Vub21pY19wYWlyZWRfZmlsZSA8LSBmaWxlLnBhdGgoZmlsZXNfZGlyLCAiZ2Vub21pY19hc3NheXNfbWF0Y2hlZF90aW1lX3BvaW50cy50c3YiKQp0bWJfdmFmX2ZpbGUgPC0gZmlsZS5wYXRoKHJlc3VsdHNfZGlyLCAidG1iX3ZhZl9nZW5vbWljLnRzdiIpCnBhbGV0dGVfZmlsZSA8LSBmaWxlLnBhdGgocm9vdF9kaXIsICJmaWd1cmVzIiwgInBhbGV0dGVzIiwgIm9uY29wcmludF9jb2xvcl9wYWxldHRlLnRzdiIpCgojIEZpbGUgcGF0aCB0byBwbG90IGRpcmVjdG9yeQpwbG90c19kaXIgPC0KICBmaWxlLnBhdGgoYW5hbHlzaXNfZGlyLCAicGxvdHMiKQppZiAoIWRpci5leGlzdHMocGxvdHNfZGlyKSkgewogIGRpci5jcmVhdGUocGxvdHNfZGlyKQp9Cgpzb3VyY2UocGFzdGUwKHJvb3RfZGlyLCAiL2ZpZ3VyZXMvc2NyaXB0cy90aGVtZS5SIikpCnNvdXJjZShwYXN0ZTAoYW5hbHlzaXNfZGlyLCAiL3V0aWwvZnVuY3Rpb24tY3JlYXRlLWJhcnBsb3QuUiIpKQpgYGAKCiMgUmVhZCBpbiBkYXRhIGFuZCBwcm9jZXNzCgpgYGB7ciBsb2FkLXByb2Nlc3MtaW5wdXRzfQpwYnRhX2RmIDwtIHJlYWRyOjpyZWFkX3RzdihwYnRhX2ZpbGUsIGd1ZXNzX21heCA9IDEwMDAwMCwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkgJT4lIAogIHNlbGVjdChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lELCBjZ19tdWx0aXBsZSwgY2dfaWQsIGNnR0ZBQywgdHVtb3JfZGVzY3JpcHRvcikKCnRtYl92YWZfZGYgPC0gcmVhZHI6OnJlYWRfdHN2KHRtYl92YWZfZmlsZSwgZ3Vlc3NfbWF4ID0gMTAwMDAwLCBzaG93X2NvbF90eXBlcyA9IEZBTFNFKSAlPiUgCiAgZmlsdGVyKCF0bWIgPj0gMTApICU+JSAKICBzZWxlY3QoS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgVmFyaWFudF9DbGFzc2lmaWNhdGlvbiwgZ2VuZV9wcm90ZWluLCBtdXRhdGlvbl9jb3VudCwJcmVnaW9uX3NpemUsIHRtYiwgVkFGKQoKZ2Vub21pY19wYWlyZWRfZGYgPC0gcmVhZHI6OnJlYWRfdHN2KGdlbm9taWNfcGFpcmVkX2ZpbGUsIGd1ZXNzX21heCA9IDEwMDAwMCwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkgJT4lCiAgbGVmdF9qb2luKHBidGFfZGYsIGJ5ID0gYygiS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCIpKSAlPiUgCiAgbGVmdF9qb2luKHRtYl92YWZfZGYsIGJ5ID0gYygiS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCIpKSAlPiUKICBmaWx0ZXIoIWlzLm5hKHRtYikpCgojIEF0dGVudGlvbiBhcyBzb21lIGJzIHNwZWNpbWVuIG1pZ2h0IG5vdCBoYXZlIFRNQiEKIyBJZiB0aGF0IGhhcHBlbnMsIHdlIHdpbGwgZW5kIHVwIHdpdGggc2FtcGxlcyBsYWNraW5nIHRpbWVwb2ludHMuCgojIFdoaWNoIHBhdGllbnQgc2FtcGxlcyBkb24ndCBoYXZlIFRNQj8KIyBnZW5vbWljX3BhaXJlZF9kZiAlPiUgCiMgIGZpbHRlcihpcy5uYSh0bWIpKSAlPiUgCiMgIHVuaXF1ZSgpICU+JSAKIyAgcmVndWxhcnRhYmxlKCkgJT4lCiMgIGZvbnRzaXplKHNpemUgPSAxMiwgcGFydCA9ICJhbGwiKQoKZGVzY3JpcHRvcnNfZGYgPC0gZ2Vub21pY19wYWlyZWRfZGYgJT4lCiAgZ3JvdXBfYnkoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkgJT4lCiAgc3VtbWFyaXplKGRlc2NyaXB0b3JzID0gcGFzdGUoc29ydCh0dW1vcl9kZXNjcmlwdG9yKSwgY29sbGFwc2UgPSAiLCAiKSwpIAoKIyBWZWN0b3IgdG8gb3JkZXIgdGltZXBvaW50cwp0aW1lcG9pbnRzIDwtIGMoIkRpYWdub3NpcyIsICJQcm9ncmVzc2l2ZSIsICJSZWN1cnJlbmNlIiwgIkRlY2Vhc2VkIiwgIlNlY29uZCBNYWxpZ25hbmN5IiwgIlVuYXZhaWxhYmxlIikKCiMgQ3JlYXRlIGRmIHRvIHVzZSBmb3IgcGxvdHMKZGZfcGxvdCA8LSBnZW5vbWljX3BhaXJlZF9kZiAlPiUgCiAgbGVmdF9qb2luKGRlc2NyaXB0b3JzX2RmLCBieSA9IGMoIktpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQiLCAiZGVzY3JpcHRvcnMiKSkgJT4lIAogIG11dGF0ZSh0ZF9jZ0dGQUMgPSBjYXNlX3doZW4oZ3JlcGwoIkRlY2Vhc2VkIiwgdHVtb3JfZGVzY3JpcHRvcikgfiAieERlY2Vhc2VkIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIFRSVUUgfiB0dW1vcl9kZXNjcmlwdG9yKSwKICAgICAgICAgbWF0Y2hfaWQgPSBwYXN0ZSh0dW1vcl9kZXNjcmlwdG9yLCBLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCBzZXAgPSAiXyIpLAogICAgICAgICB0dW1vcl9kZXNjcmlwdG9yID0gZmFjdG9yKHR1bW9yX2Rlc2NyaXB0b3IpLAogICAgICAgICB0dW1vcl9kZXNjcmlwdG9yID0gZmN0X3JlbGV2ZWwodHVtb3JfZGVzY3JpcHRvciwgdGltZXBvaW50cykpCgojIExldCdzIGNvdW50IG51bWJlciBvZiBzYW1wbGVzIApjb3VudF9kZiA8LSBkZl9wbG90ICU+JSAKICBncm91cF9ieSh0dW1vcl9kZXNjcmlwdG9yLCBjZ19pZCwgbWF0Y2hfaWQsIFZhcmlhbnRfQ2xhc3NpZmljYXRpb24pICU+JSAKICBkcGx5cjo6Y291bnQoY2dfaWQpIAoKYGBgIAoKIyBEZWZpbmUgcGFyYW1ldGVycyBmb3IgcGxvdHMKCmBgYHtyIGRlZmluZS1wYXJhbWV0ZXJzLWZvci1wbG90c30KIyBSZWFkIGNvbG9yIHBhbGV0dGUKcGFsZXR0ZV9kZiA8LSByZWFkcjo6cmVhZF90c3YocGFsZXR0ZV9maWxlLCBndWVzc19tYXggPSAxMDAwMDAsIHNob3dfY29sX3R5cGVzID0gRkFMU0UpIAoKIyBEZWZpbmUgYW5kIG9yZGVyIHBhbGV0dGUKcGFsZXR0ZSA8LSBwYWxldHRlX2RmJGhleF9jb2RlcwpuYW1lcyhwYWxldHRlKSA8LSBwYWxldHRlX2RmJGNvbG9yX25hbWVzCmBgYAoKIyBBbHRlcmF0aW9ucyBwZXIgdGltZXBvaW50CgpgYGB7ciBwbG90LXRpbWVwb2ludCwgZmlnLndpZHRoID0gNiwgZmlnLmhlaWdodCA9IDYsIGZpZy5mdWxsd2lkdGggPSBUUlVFfQojIERlZmluZSBwYXJhbWV0ZXJzIGZvciBmdW5jdGlvbgp4X3ZhbHVlIDwtIGNvdW50X2RmJHR1bW9yX2Rlc2NyaXB0b3IKdGl0bGUgPC0gcGFzdGUoIlZhcmlhbnQgdHlwZXMgaW4gUEJUQSBjb2hvcnQiLCBzZXAgPSAiICIpCgojIFJ1biBmdW5jdGlvbgpmbmFtZSA8LSBwYXN0ZTAocGxvdHNfZGlyLCAiLyIsICJBbHRlcmF0aW9uX3R5cGVfdGltZXBvaW50c19iYXJwbG90cy5wZGYiKQpwcmludChmbmFtZSkKcCA8LSBjcmVhdGVfc3RhY2tlZF9iYXJwbG90X3ZhcmlhbnQoY291bnRfZGYgPSBjb3VudF9kZiwgeCA9IHhfdmFsdWUsIHBhbGV0dGUgPSBwYWxldHRlLCB0aXRsZSA9IHRpdGxlKQpwZGYoZmlsZSA9IGZuYW1lLCB3aWR0aCA9IDYsIGhlaWdodCA9IDYpCnByaW50KHApCmRldi5vZmYoKQpgYGAKCiMgQWx0ZXJhdGlvbnMgcGVyIHRpbWVwb2ludCBpbiBlYWNoIGNhbmNlciB0eXBlCgpgYGB7ciBwbG90LWNnLWlkLCBmaWcud2lkdGggPSAzMCwgZmlnLmhlaWdodCA9IDE4LCBmaWcuZnVsbHdpZHRoID0gVFJVRX0KY2dfaWRfaWQgPC0gYXMuY2hhcmFjdGVyKHVuaXF1ZShjb3VudF9kZiRjZ19pZCkpIApjZ19pZF9pZCA8LSBzb3J0KGNnX2lkX2lkLCBkZWNyZWFzaW5nID0gRkFMU0UpCmNnX2lkX2lkCgojIERlZmluZSBwYXJhbWV0ZXJzIGZvciBmdW5jdGlvbgp4X3ZhbHVlIDwtIGNvdW50X2RmJG1hdGNoX2lkCnRpdGxlIDwtIHBhc3RlKCJWYXJpYW50IHR5cGVzIGluIFBCVEEgY29ob3J0IGFjcm9zcyBjYW5jZXIgZ3JvdXBzIiwgc2VwID0gIiAiKQoKIyBSdW4gZnVuY3Rpb24KZm5hbWUgPC0gcGFzdGUwKHBsb3RzX2RpciwgIi8iLCAiQWx0ZXJhdGlvbl90eXBlX3RpbWVwb2ludHNfY2dfaWRfYmFycGxvdHMucGRmIikKcHJpbnQoZm5hbWUpCnAgPC0gY3JlYXRlX3N0YWNrZWRfYmFycGxvdF92YXJpYW50X2NnX2lkKGNvdW50X2RmID0gY291bnRfZGYsIHggPSB4X3ZhbHVlLCBwYWxldHRlID0gcGFsZXR0ZSwgdGl0bGUgPSB0aXRsZSkKcGRmKGZpbGUgPSBmbmFtZSwgd2lkdGggPSAzMCwgaGVpZ2h0ID0gMTgpCnByaW50KHApCmRldi5vZmYoKQoKYGBgCgojIEFsdGVyYXRpb25zIHBlciB0aW1lcG9pbnQgaW4gZWFjaCBjYW5jZXIgdHlwZSBkZWZpbmVkIGJ5IGNnR0ZBQwoKYGBge3IgcGxvdC1jZ0dGQUMtbi1pbmRpdmlkdWFsLXBsb3RzLCBmaWcud2lkdGggPSA1LCBmaWcuaGVpZ2h0ID0gNSwgZmlnLmZ1bGx3aWR0aCA9IFRSVUV9CmRmX3Bsb3RfY2dHRkFDIDwtIGRmX3Bsb3QgJT4lCiAgZ3JvdXBfYnkodHVtb3JfZGVzY3JpcHRvciwgY2dHRkFDLCBtYXRjaF9pZCwgVmFyaWFudF9DbGFzc2lmaWNhdGlvbikgJT4lIAogIGRwbHlyOjpjb3VudChjZ0dGQUMpCgpjZ0dGQUNfaWQgPC0gYXMuY2hhcmFjdGVyKHVuaXF1ZShkZl9wbG90X2NnR0ZBQyRjZ0dGQUMpKSAKY2dHRkFDX2lkIDwtIHNvcnQoY2dHRkFDX2lkLCBkZWNyZWFzaW5nID0gRkFMU0UpCmNnR0ZBQ19pZAoKIyBEZWZpbmUgcGFyYW1ldGVycyBmb3IgZnVuY3Rpb24KeF92YWx1ZSA8LSBkZl9wbG90X2NnR0ZBQyRtYXRjaF9pZAp0aXRsZSA8LSBwYXN0ZSgiVmFyaWFudCB0eXBlcyBpbiBQQlRBIGNvaG9ydCBhY3Jvc3MgY2dHRkFDIiwgc2VwID0gIiAiKQoKIyBSdW4gZnVuY3Rpb24KcCA8LSBjcmVhdGVfc3RhY2tlZF9iYXJwbG90X3ZhcmlhbnRfY2dHRkFDKGNvdW50X2RmID0gZGZfcGxvdF9jZ0dGQUMsIHggPSB4X3ZhbHVlLCBwYWxldHRlID0gcGFsZXR0ZSwgdGl0bGUgPSB0aXRsZSkKCmBgYAoKIyBBbHRlcmF0aW9ucyBwZXIgdGltZXBvaW50IGluIGVhY2ggY2FuY2VyIHR5cGUgYW5kIHRpbWVwb2ludCBtb2RlbAoKYGBge3IgcGxvdC10aW1lcG9pbnQtbW9kZWwsIGZpZy53aWR0aCA9IDE1LCBmaWcuaGVpZ2h0ID0gMTIsIGZpZy5mdWxsd2lkdGggPSBUUlVFfQp0bV9kZl9wbG90IDwtIGRmX3Bsb3QgJT4lCiAgZmlsdGVyKCFpcy5uYSh0aW1lcG9pbnRzX21vZGVscykpICU+JSAKICBncm91cF9ieSh0dW1vcl9kZXNjcmlwdG9yLCBjZ19pZCwgdGltZXBvaW50c19tb2RlbHMsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIFZhcmlhbnRfQ2xhc3NpZmljYXRpb24pICU+JSAKICBkcGx5cjo6Y291bnQodGltZXBvaW50X2NnR0ZBQ19uKQoKdG0gPC0gYXMuY2hhcmFjdGVyKHVuaXF1ZSh0bV9kZl9wbG90JHRpbWVwb2ludHNfbW9kZWxzKSkKdG0gPC0gc29ydCh0bSwgZGVjcmVhc2luZyA9IEZBTFNFKQp0bQoKIyBMb29wIHRocm91Z2ggdmFyaWFibGUKZm9yIChpIGluIHNlcV9hbG9uZyh0bSkpewogIHByaW50KGkpCiAgZGZfc3ViIDwtIHRtX2RmX3Bsb3QgJT4lCiAgICAgIGZpbHRlcih0aW1lcG9pbnRzX21vZGVscyA9PSB0bVtpXSkKICAKICAgIyBEZWZpbmUgcGFyYW1ldGVycyBmb3IgZnVuY3Rpb24KICB4X3ZhbHVlIDwtIGRmX3N1YiRtYXRjaF9pZAogIHRpdGxlIDwtIHBhc3RlKHRtW2ldKQogIAogICMgUnVuIGZ1bmN0aW9uCiAgcCA8LSBjcmVhdGVfc3RhY2tlZF9iYXJwbG90X3ZhcmlhbnRfY2dfaWQoY291bnRfZGYgPSBkZl9zdWIsIHggPSB4X3ZhbHVlLCBwYWxldHRlID0gcGFsZXR0ZSwgdGl0bGUgPSB0aXRsZSkKICAKfQpgYGAKCgpgYGB7ciBlY2hvPVRSVUV9CnNlc3Npb25JbmZvKCkKYGBgCg==